Boundary conditions: The path integral approach 
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Abstract. The path integral approach to quantum mechanics requires a substantial 
generalisation to describe the dynamics of systems confined to bounded domains. Non- 
local boundary conditions can be introduced in Feynman's approach by means of boundary 
amplitude distributions and complex phases to describe the quantum dynamics in terms of 
the classical trajectories. The different prescriptions involve only trajectories reaching the 
boundary and correspond to different choices of boundary conditions of selfadjoint extensions 
of the Hamiltonian. One dimensional particle dynamics is analysed in detail. 



1. Introduction 

The original approach to quantum mechanics due to Heisenberg and Schrodinger is based on 
the Hamiltonian formalism. Dirac and later Feynman introduced a Lagrangian approach. 
The Feynman's path integral method was shown later to be equivalent to the traditional 
formulation for most of quantum mechanical systems [1, 2]. In field theory the method presents 
some advantages over Hamiltonian quantisation. The Lagrange formalism preserves relativistic 
covariance which makes the Feynman method very convenient to achieve the renormalization of 
field theories both in perturbative and non-perturbative approaches. In the Euclidean version 
[3j the functional integral formulation is a crucial ingredient for non-perturbative approaches to 
field theory and critical phenomena [3]. The equivalence of the functional integral formulation 
with the Hamiltonian pictures also holds for constrained systems like gauge theories and string 
theories. 

However, for systems constrained to bounded domains Feynman's approach becomes more 
intricate than the Hamiltonian operator approach and require a sophisticated implementation for 
generic types of boundary conditions compatible with Hamiltonian approach [5j. The analysis 
of this problem is the main goal of this paper. On the other hand the physics of quantum 
systems constrained to bounded domains is becoming very relevant not only for applications in 
condensed matter (quantum Hall efi'ect, graphene, quantum wires, quantum dots, etc) but also 
in fundamental physics in fields like quantum gravity and string theory. In all these phenomena 
the role of boundary conditions is very important to describe a variety of new physical eff'ects like 
anomalies [El El [8] topology change ^ quantum holography [TOl HH 12] , quantum gravity and 
AdS/CFT correspondence [13j . To some extent the relevance of boundaries in the description of 
fundamental physical phenomena, has promoted the role of boundary phenomena from academic 
and phenomelogical simplifications of complex physical systems to a higher status connected 
with very basic fundamental principles. The application of the path integral approach to those 
problems requires a radical generalisation of Feynman's formalism. 



2. Boundary conditions in quantum mechanics 



Unitarity is the basic quantum principle which governs the dynamics of quantum systems 
confined in bounded domains. The analytical implementation of this condition in the selfadjoint 
character of the Hamiltonian operator encodes all the quantum subtleties associated to the 
unitarity principle and the dynamical behaviour at the boundary. The selfadjointness condition 
of the Hamiltonian imposes severe restrictions to the boundary conditions which are compatible 
with quantum mechanics. 

The existence of a boundary generically enhances the genuine quantum aspects of the system. 
Well known examples of this enhancement are Young's two slits experiments and the Aharonov- 
Bohm effect, which pointed out the relevance of boundary conditions in the quantum theory. 
Another examples of quantum physical phenomena which are intimately related to boundary 
conditions are the Casimir effect [14:\ [T5] the role of edge states |16j and the quantisation of 
conductivity [T71 [18] in the quantum Hall effect. 

The quantum role of boundary conditions is very important for the behaviour of low energy 
levels. High energy levels are quite independent of boundary effects. Indeed, boundary effects 
play no role in the ultraviolet regime, whereas they are crucial for the infrared |19j . 

Let us consider a point-like particle moving on a bounded domain Q of M."' with smooth 
oriented boundary dO,. The Hamiltonian of the system can be constructed from the scalar 
Laplacian 



This operator is symmetric on the domain of smooth wave functions with compact support on fi. 
However, the restriction of the Laplacian operator to this domain is not selfadjoint, but it can 
be extended to larger domains of wave functions where it becomes selfadjoint. The extension, 
however, not unique. 

The classification of selfadjoint extensions of the Hamiltonian can be characterised in terms of 
unitary operators between defect subspaces in the classical theory due to von Neumann [201 [2T]. 
However, there is a more useful characterisation of these selfadjoint extensions in terms of 
constraint conditions on the boundary values of physical wave functions [22]. In this framework 
the set of self-adjoint extensions of the Hamiltonian is in one-to-one correspondence with the 
group of unitary operators of the Hilbert space of wave functions of the boundary 517 

which are square integrable with respect to the Riemaniann measure induced from the Euclidean 
metric of M."'. 

Thus, any unitary operator U of the Hilbert space L^{dQ) defines a selfadjoint quantum 
Hamiltonian A^. And conversely, any selfadjoint extension of H is associated to one unitary 
operator U of this type [22] . The domain of the selfadjoint Hamiltonian governed by U is defined 
by the wave functions which satisfy the boundary condition 



where (p = iplg^ is the boundary value of the wave function ^ and its oriented normal derivative 
at the boundary dil. Through the above characterisation, the set of self-adjoint extensions of 
the Hamiltonian inherits the group structure of the group of unitary operators. For spaces of 
dimension higher than one the group of boundary conditions is an infinite dimensional group. 

Since we shall focus on one-dimensional spaces two domains are of special interest: the 
compact interval = [0,1] and the half line = [O.oo). Although the second domain is non- 
compact it will be useful to illustrate the main problems due to its simplicity. In the first case 
the set of boundary conditions is U{2) whereas for the half line it is U{1). 





i=l 



ip - iip = U{ip + iip), 



(2) 



3. Boundary conditions in path integrals 



In Feynman's approach to quantum mechanics the dynamics is governed by an action principle 
similar to that which governs the classical dynamics. Whereas the classical dynamics is given, 
according to the variational action principle, by stationary trajectories of the classical action, 
the quantum dynamics is automatically implemented in the path integral formalism by the 
weight that the classical action provides for particle trajectories. In the Euclidean time T = —it 
formalism the evolution propagator for the free particle is given by the path integral 

K,{x,y) = e-^^(x,y) = / [5x{t)] e'^o -W^-^*. (3) 

J x{T) = y 

However, for particles evolving in a bounded domain the variational problem is not uniquely 
defined in classical mechanics. In order to have a deterministic evolution it becomes necessary to 
specify the evolution of the particles after reaching the boundary. A similar ambiguity problem 
arises in quantum mechanics. 

The boundary imposes more severe constraints on the classical dynamics than in the quantum 
one. This is due to the point-like nature of the particle which requires that after reaching the 
boundary the individual particle has to reemerge back either at the same point or at a different 
point of the boundary. In pure classical mechanics the only freedom concerns where it reemerges 
back and the momentum it reemerges with. The emergence of the particle at a different point 
grants the possibility of folding and gluing the boundary of the domain giving rise to non- 
trivial topologies. In summary, the classical boundary conditions are given by two maps [23j: 
an isometry of the boundary a : d^l — > dQ and a positive density function p : d^l M"*" 
which specify the change of position and normal component of momentum of the trajectory 
of the particle upon reaching the boundary. The isometry a encodes the possible geometry 
and topology generated by the folding of the boundary and the function p is associated to 
the reflectivity (transparency or stickiness) properties of the boundary. However, the quantum 
boundary conditions have a larger set of possibilities described by a unitary group. In order to 
have a path integral description of all boundary conditions we need to incorporate some random 
behaviour for the trajectories reaching the boundary and complex phases for those trajectories. 
This is possible because the wave functions are complex and the evolution operator involves 
complex amplitudes. Although, is this way we are able to describe any type of unitary evolution 
in the bounded domain the method goes far beyond Feynman's pure action approach. 

The prescription is quite involved and proceeds by considering instead of the Euclidean time 
evolution propagator K^, the resolvent operator Cz of the Hamiltonian 

fOO irp 

C,{x,y) = {zl +Hr\x,y) = —^-^^ K^{x,y). (4) 

The Euclidean time propagator can be recovered from the resolvent by means of the following 
countour integral 

K^{x,y) = ^.j Cz{x,yy^dz (5) 

along a contour which encloses the spectrum of H on the real axis. 

Boundary conditions can be easily implemented into the resolvent, whereas as we shall see 
that the implementation in the Euclidean time propagator is much harder. Let us consider a 
fixed boundary condition, e.g. the Neumann boundary conditions Uq = 1 , and consider the 
corresponding Hamiltonian as a background selfadjoint operator. The selfadjoint extension 
of H defined on the domain 

i{l + U)^ = (I - U)<f (6) 



by the unitary operator U has a resolvent given by Krein's formula |24j 

Cf(x,y) = CO(x,y)- / dw[ dw' Cl{x,w)R'i {wM)Cl{w' ,y), (7) 
Jan Jan 

where is the operator of L'^{dQ) defined by 

R'^ = ((I - - z(/ + U)r\l - U). (8) 

A similar formula could be obtained choosing another boundary condition as background 
boundary condition instead of Neumann's condition. 

The inverse transform permits to recover a formula for the propagator kernel of the type 

K^{x,y)=Kl{x,y)-^ I dze^^ f dwf dw' Cl{x,w)R'i {w,w')Cl{w' ,y). (9) 

It is easy to rewrite K^{x,y) as a path integral as in ([3]) restricting the trajectories to the 
interior of the domain and counting twice the trajectories hitting the boundary dO,. However, 
in general, the kernel K^{x,y) cannot be rewritten in terms of a path integral. Only for a 
few boundary conditions the reduction can be achieved, but for generic boundary conditions 
the kernel K^{x,y) has to be considered as a genuine boundary condition kernel containing 
information about the boundary jumps amplitudes and phases associated to the different 
trajectories hitting the boundary. The complex structure of this kernel reduces the utility of 
the path integral approach and points out the behaviour of the boundary as a genuine quantum 
device. This behaviour can be explicitly pointed out by noticing that under certain boundary 
conditions the quantum evolution of a narrow wave packet is scattered backward by the boundary 
as a quite widespread wave packet emerging from all points of the boundary. 

However, there are cases where this kernel adopts a simple form and the path integral 
approach can be formulated in very explicit way. In particular, for Dirichlet boundary conditions 
C/ = -I , 

R^i = {cir' (10) 

and 

C^{x,y)=Cl{x,y)- f dwf dw' Cl{x,w){C^,)-\w,w')C^Aw\y), (11) 
Ja^ Jan 

which leads to a propagator kernel given by the path integral ^ but restricted to paths which 
do not reach the boundary dil. 

Further examples can be explicitly analysed for one-dimensional domains. 



4. Particle in a half line 

In order to illustrate the method let us consider for simplicity the unbounded domain = [0, cxd]. 
The set of selfadjoint extensions is parametrised by U{1), i.e the unitary operators of the form 
U = e*". The corresponding boundary conditions are the mixed conditions 

(^(0) = V(0) = - tan -99(0). (12) 

The resolvent of the corresponding selfadjoint extension can be obtained from the Krein formula 
d?]) giving rise to 



Cnx,y) 



1 



-V2z\x-y\ _|_ 



1 2z 



tanf + 



2ztan| 



'2z 



(13) 



which in the two extreme cases a = (Neumann) 



CO(x, y) = -^e-^A-y\ + ^^^-V2-z(\A+\y\) (14) 



and a = § (Dirichlet) 

Cll\x, y) = J—^-^A^-y\ _ J_^-^-{\A+\y\) (15) 

is reduced to very simple formulas. Notice that because 

iif = V2i (16) 

the Dirichlet resolvent Cj''^ can be easily obtained from the Neumann resolvent from Krein's 
formula ([7]). 

The inverse Laplace transform gives the following formulas for the Euclidean evolution kernels 
with Neumann boundary conditions [25t i26j 



and Dirichlet boundary conditions 

^ ^/2^ ^2^ 

which can be interpreted in terms of the path integral as a sum over paths which never reach 
the boundary [25^ I26j. whereas the Neumann kernel corresponds to double-counting the paths 
which hit the boundary 0. Notice that the terms depending on x + y, which appear in ()17p 
(lisp (I14p and ()15p . break translation invariance due to the existence of finite boundaries of the 
interval. 



5. Particle in an interval. 

If the particle is confined in an interval [0, 1] the set of selfadjoint extensions is parametrised by 
U{2). Although the general theory can be developed on the same basis, we shall consider only 
few cases where the path integral approach is simplified, 
i) Neumann boundary conditions U = I , (f = 0. 

The resolvent with Neumann boundary conditions at both ends of the interval is (see e.g. 

[271 [2H]) 

I _ Q-2V2Z I _ g-2^/2i ' ~^ I _ I _ Q-2V2Z 

The Euclidean time propagator kernel 



1 



2V2z 



-(x+y+m)'^/2T 



^ In the Euclidean approach the restrictions imposed on the paths for Neumann and Dirichlet boundary conditions 
are interchanged with respect to the boundary conditions in the Minkowski approach, where the Dirichlet problem 
is associated to an infinite totally refiecting wall. 



corresponds to a path integral where the trajectories which hit the boundary are double weighted 

as in the half line case. 

ii) Dirichlct boundary conditions U = —I , ip{0) = '/'(I) = 

The resolvent with Dirichlet boundary conditions at both ends of the interval is 

I _ g-2^2z I _ Q-2V2Z I _ Q-2V2Z I _ Q-2V2Z 

The associated propagator kernel 



2V2z 



00 00 
12, 



-{x-y+my/2T _ 



-ix+y+m)^ /2T 



corresponds to a path integral where the trajectories reaching the boundary have been removed. 

^0 r 



iii) Periodic boundary conditions U = ai 

The resolvent 



and the propagator kernel 



1 



2V2z 



-V2^\x-y\ _^ 6 
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2z{l-x+y) 



2z(l+x-y) 



\ _ Q-\/2z _ Q—\/2z 



(19) 



-. 00 



indicate that the corresponding path integral is performed over periodic trajectories which cross 
from one end of the interval to the other, 
iv) Pseudo-periodic boundary conditions. 

The selfadjoint Hamiltonian corresponding to the unitary operator 



U 



e"' 
e*^ 



is defined in the domain satisfying the boundary conditions 

(^(1) = e-(^(0) (^'(l) = eV(0). 

The resolvent 



and propagator 



Je{y-x) 



2V2z 



-V2i\x-y\ _(_ 6 



2z(l-x+y) 



1 - e-^- 



-V2z(l-y+x) 

- e-'' + ^— e'' 



]^ Q—V2z+ie 



(20) 



(21) 



(22) 



^ 00 

(x v) = -{x-y+nf /2T+ie{n+x-y) 



kernels correspond to a path integral over periodic trajectories which are weighted with and 
additional phase factor e**^ for everytime they cross the (periodic) boundary from left to right 
and e~*^ if they cross it in opposite direction. 



The method of images also permits to use unconstrained path integrals to describe systems 
with non-trivial boundary conditions [29j . However, in the case higher dimensions the method is 
not useful in the presence of non symmetric boundaries and the path integral cannot be defined 
by a simple prescription as in the Feynman original formulation. 

For some boundary conditions of one dimensional systems it is possible to use another method 
based on the path integral approach with simpler boundary conditions but with an additional 
singular potential in the Euclidean action [30J. Two interesting cases are the following. 

v) Quasi-periodic boundary conditions. 

The selfadjoint extension corresponding to the unitary operator 

jj _ f cos a sin a 
" ysina —cos a 

is defined in the domain of quasiper iodic wave functions 

ip{l) = tan - vp(0); ^' {I) = cotan- ^^'(0). (24) 

In this case the resolvent, Euclidean time propagators and path integral can be identified with 
those of periodic boundary conditions with an extra potential term in the action corresponding 
to a 5' {x) singular interaction. 

vi) A circle with a defect point. 

The Hamiltonian with boundary conditions corresponding to the unitary matrix 

U = ^('^ .^1 (25) 
2-ia\2 laj ^ ' 

can be thought as equivalent to the Hamiltonian of a delta function a5{x) potential acting on 
a particle moving a circle with the usual periodic boundary conditions U = ai \31\ |32] The 
corresponding formulae for the propagator kernel can be also written in a compact way \32\ I33j. 

However, the method is only restricted to similar cases and for generic boundary conditions a 
closed form expression is not available. In higher dimensions the number of boundary conditions 
for which the path integral method is useful to describe the quantum evolutions is even more 
limited. 

6. Conclusions 

In summary, it is possible to generalise the Feynman approach to describe the dynamics of 
quantum systems constrained to bounded domains. The boundary itself has to be considered 
from this point of view a genuine quantum device and transitions amplitudes are required to 
implement the effect of the boundary into the path integral approach. For some boundary 
conditions the modification of the path integral formula includes a phase factor or a boundary 
weight for the trajectories which reach the boundary. However, the method becomes not useful 
for generic boundary conditions because the prescription becomes very intricate. This fact is a 
consequence of the enhancement of genuine quantum effects by the presence of the boundary. 



(23) 
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